NCBR-227

Cox Regression on GATA2 data

Author

Subrata Paul

Data

Kaplan-Meier Analysis

Mixed Effect Cox Regression

Cox proportional hazard model is commonly used in survival analysis that models the hazard function which denotes the instantaneous rate of the occurrence of an event at a time given that the subject survived upto that point. Including random effect in the Cox’s proportional hazard model enables us to incorporate within group homogeneity in outcomes. A mixed effect model in general multivariate analysis is used to control for relatedness in the samples. I used the kinship matrix (kinship2 R package) to correct for relatedness as a random effect.

  Surv(ageonset,censor)
Predictors HR CI p
mutation Class [C term] 1.79 0.53 – 6.05 0.346
mutation Class [Null] 3.36 1.67 – 6.77 0.001
mutation Class
[Truncation]
5.70 2.42 – 13.42 <0.001
Gender [2] 1.06 0.58 – 1.92 0.852
N Indiv.ID 128
Observations 128

Enhancer is selected as the baseline mutation class. All the mutation other than C-term increases risk compared to enhancer mutation. An individual with truncation mutation has five times hazard than an individual with enhancer mutation. To be able to compare each mutation class, I used emmeans package. Note that, results in the following table are averaged over the levels of: Gender. P value adjustment was made using tukey method for comparing a family of 5 estimates.

Most important assumption for this model is the “proportional hazard” which was not violated.

               chisq df    p
mutation.Class  4.65  3 0.20
Gender          0.87  1 0.35
GLOBAL          4.86  4 0.30

Sample size for testing difference between missense group is 81. R398 is selected as the baseline missense group. Missense group R396 (strong evidence p < 0.01) and T354 (moderate evidence p< 0.05) increases risk compared to R398 missense mutation. An individual with R396 missense mutation on an average has two and half times higher hazard than an individual with R398 mutation.

  Surv(ageonset,censor)
Predictors HR CI p
Missense group [R361] 1.63 0.79 – 3.35 0.186
Missense group [R396] 2.67 1.36 – 5.26 0.004
Missense group [T354] 2.16 1.09 – 4.31 0.028
Gender [2] 1.17 0.69 – 1.99 0.549
N Indiv.ID 81
Observations 81

Packages and Citations

R version 4.4.2 (2024-10-31)
Platform: aarch64-apple-darwin20
Running under: macOS Ventura 13.7.4

Matrix products: default
BLAS:   /Library/Frameworks/R.framework/Versions/4.4-arm64/Resources/lib/libRblas.0.dylib 
LAPACK: /Library/Frameworks/R.framework/Versions/4.4-arm64/Resources/lib/libRlapack.dylib;  LAPACK version 3.12.0

locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8

time zone: America/New_York
tzcode source: internal

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
 [1] kableExtra_1.4.0 kinship2_1.9.6.1 quadprog_1.5-8   Matrix_1.7-1    
 [5] xtable_1.8-4     coxme_2.2-22     bdsmatrix_1.3-7  survminer_0.5.0 
 [9] ggpubr_0.6.0     DT_0.33          lubridate_1.9.3  forcats_1.0.0   
[13] stringr_1.5.1    dplyr_1.1.4      purrr_1.0.4      readr_2.1.5     
[17] tidyr_1.3.1      tibble_3.2.1     ggplot2_3.5.1    tidyverse_2.0.0 
[21] survival_3.7-0   readxl_1.4.3    

loaded via a namespace (and not attached):
 [1] gridExtra_2.3      sandwich_3.1-1     rlang_1.1.5        magrittr_2.0.3    
 [5] multcomp_1.4-26    snakecase_0.11.1   compiler_4.4.2     systemfonts_1.1.0 
 [9] vctrs_0.6.5        pkgconfig_2.0.3    fastmap_1.2.0      backports_1.5.0   
[13] labeling_0.4.3     effectsize_0.8.9   KMsurv_0.1-5       rmarkdown_2.29    
[17] markdown_1.13      tzdb_0.4.0         nloptr_2.1.1       xfun_0.50         
[21] cachem_1.1.0       jsonlite_1.8.9     sjmisc_2.8.10      ggeffects_1.7.2   
[25] broom_1.0.7        R6_2.5.1           bslib_0.9.0        stringi_1.8.4     
[29] boot_1.3-31        car_3.1-3          jquerylib_0.1.4    cellranger_1.1.0  
[33] estimability_1.5.1 Rcpp_1.0.14        knitr_1.49         zoo_1.8-12        
[37] parameters_0.23.0  splines_4.4.2      timechange_0.3.0   tidyselect_1.2.1  
[41] rstudioapi_0.17.1  abind_1.4-8        yaml_2.3.10        codetools_0.2-20  
[45] ggtext_0.1.2       sjlabelled_1.2.0   lattice_0.22-6     withr_3.0.2       
[49] bayestestR_0.15.0  evaluate_1.0.3     xml2_1.3.6         survMisc_0.5.6    
[53] pillar_1.10.1      carData_3.0-5      insight_0.20.5     generics_0.1.3    
[57] hms_1.1.3          munsell_0.5.1      commonmark_1.9.2   scales_1.3.0      
[61] minqa_1.2.8        glue_1.8.0         sjPlot_2.8.16      emmeans_1.10.5    
[65] tools_4.4.2        data.table_1.16.4  lme4_1.1-35.5      ggsignif_0.6.4    
[69] mvtnorm_1.3-2      grid_4.4.2         crosstalk_1.2.1    datawizard_0.13.0 
[73] colorspace_2.1-1   nlme_3.1-166       performance_0.12.4 Formula_1.2-5     
[77] cli_3.6.3          km.ci_0.5-6        viridisLite_0.4.2  svglite_2.1.3     
[81] sjstats_0.19.0     gtable_0.3.6       rstatix_0.7.2      sass_0.4.9        
[85] digest_0.6.37      TH.data_1.1-2      htmlwidgets_1.6.4  farver_2.1.2      
[89] htmltools_0.5.8.1  lifecycle_1.0.4    MASS_7.3-61        gridtext_0.1.5    

Citation for R itself

NULL

Citation for survival package

To cite package 'survival' in publications use:

  Therneau T (2024). _A Package for Survival Analysis in R_. R package
  version 3.7-0, <https://CRAN.R-project.org/package=survival>.

  Terry M. Therneau, Patricia M. Grambsch (2000). _Modeling Survival
  Data: Extending the Cox Model_. Springer, New York. ISBN
  0-387-98784-3.

To see these entries in BibTeX format, use 'print(<citation>,
bibtex=TRUE)', 'toBibtex(.)', or set
'options(citation.bibtex.max=999)'.

Citation for coxme

To cite package 'coxme' in publications use:

  Therneau TM (2024). _coxme: Mixed Effects Cox Models_. R package
  version 2.2-22, <https://CRAN.R-project.org/package=coxme>.

A BibTeX entry for LaTeX users is

  @Manual{,
    title = {coxme: Mixed Effects Cox Models},
    author = {Terry M. Therneau},
    year = {2024},
    note = {R package version 2.2-22},
    url = {https://CRAN.R-project.org/package=coxme},
  }

Citation for kinship2

To cite package 'kinship2' in publications use:

  Sinnwell J, Therneau T (2024). _kinship2: Pedigree Functions_. R
  package version 1.9.6.1,
  <https://CRAN.R-project.org/package=kinship2>.

A BibTeX entry for LaTeX users is

  @Manual{,
    title = {kinship2: Pedigree Functions},
    author = {Jason Sinnwell and Terry Therneau},
    year = {2024},
    note = {R package version 1.9.6.1},
    url = {https://CRAN.R-project.org/package=kinship2},
  }